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ABSTRACT 

We present the discovery of an absorption- line redshift of z = 2.609 for GRB 090426, 
establishing the first firm lower limit to a redshift for a gamma-ray burst with an 
observed duration of < 2 s. With a rest-frame burst duration of T 90z = 0.35 s and a 
detailed examination of the peak energy of the event, we suggest that this is likely (at 
>90% confidence) a member of the short/hard phenomenological class of GRBs. From 
analysis of the optical-afterglow spectrum we find that the burst originated along a 
very low HI column density sightline, with N m < 3.2 x 10 19 cm" 2 . Our GRB 090426 
afterglow spectrum also appears to have weaker low-ionisation absorption (Si II, C 
II) than ~95% of previous afterglow spectra. Finally, we also report the discovery of 
a blue, very luminous, star-forming putative host galaxy (~ at a small angular 
offset from the location of the optical afterglow. We consider the implications of this 
unique GRB in the context of burst duration classification and our understanding of 
GRB progenitor scenarios. 

Key words: gamma-rays: bursts — galaxies: interstellar medium 



1 INTRODUCTION 

The early emergence of a two-class bifurcation in the 
high-energy properties of gamma-ray bursts (GRBs; 
iKouveliotou et al l Il993h gave rise to the supposition that 
two distinct "progenitors" could be responsible for t he lion's 
share of such events fe.g.. IZhang fc Me szaros 2004J). While 
observations directly link long-duration soft-s pectra GRBs 
(LSB s ) to the death of y oung massive stars l|Stanek et al.l 
120031 : IHiorth et al] 12003 : see IWooslev fc Blooml |2006| "for 
a review), less-strong circumstantial evidence (based on 
physical associations with more evolved galaxies) suggests 
that at least some fraction of short-duration hard-spectra 
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GRBs (SHBs) are due to older progenitors dBloom et all 
20061: iFox et al.ll2005l: iBerger et alll2005l: IHiorth et alll2005l : 



Prochaska et al . 200d: lBloom fc Prochas ka 20060. Whether 
SHBs are due to the coalescence of two neutron stars, other 
compact degenerate binaries, some combination of these 



In this paper, we will use the term "SHB" to denote this phe- 
nomenological class of gamma-ray bursts as first identified in the 
BATSE sample. This is distinct from a simple cut on duration 
or hardness, since the populations are known to overlap to vary- 
ing degrees with different instruments. Furthermore, identifica- 
tion with a class does not necessarily imply identification with 
a particular progenitor, even though most LSBs have been asso- 
ciated with massive stars and a few SHBs have been associated 
with old populations. 
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models, or something else entirely is an open question (cf. 
iLee fc Ramirez-Ruizll2007h . 

Several lines of evidence now suggest that the true pro- 
genitor diversity does not map with one-to-one correspon- 
dence to the two-class phenomenological landscapfl In par- 
ticular, there appear to be many more than just two pro- 
genitors. For example, a small fraction of SHBs probably 
originate from massive flaring activity o f extragalactic mag- 
netars (highly magnetised neutron s tars; lAbbott et al.|[200c1 ; 
iHurlevlHooi : IChapman et al.l 120091 ). Either similar magne- 
tar activity (or perhaps a flaring accretion-powered sys- 
tem) from objects in our own Galaxy may occasionally cre- 
ate L SBs as well ()Kasliwal et al . 2008; Castro -Tirad o et al.l 
2008). Classification of individual events even among the 
two well-established cosmological groups has proven ex- 
tremely difficult, not only in the overlap region of the 
duration-hardness diagram where the population distribu- 
tions merge (at around 1—2 seconds for Swift/HAT) but 
even for much longer-duration events. Indeed, some of the 
same SHBs which have been used to link this phenomeno- 
logical class to older, evolved galaxies actually have ob- 
served total durations (as measured by Tgo, the interval 
over which 90% of the burst counts are observed) of over 
100 s due to a component of extended emission (EE) that 
follows the initial spike. At least two LSBs at low red- 
shift, GRB 060505 and GRB 060614, were not accompanied 
by obse rvable supernovae despite intense follo w-up cam- 
paigns l|Fvnbo et al.l 120061 : iGehrels et ai1l2006h . and it is 
still debated whether these events group most naturally 
with short-duration events, long-duration events, or a new 
class entirely (e.g. Ijakobsson fc Fvnbo!l2007l: iLu et aD | 200S ; 
IXu et al.|[2009l ; iLevesque fc Kewlevll2007l ; IThone et al.ll200S , 
and others). Most recently, the two highest redshift GRBs 
detected to date, at z = 6.7 and z = 8.2 llGreiner et al 
20091 ; ISalvaterra et al.ll2009l ; iTanvir et al]|2009l ; IZhang et al 



2009), were observed to have rest-frame durations of < 2 s, 
yet few have argued that these events did not arise from 
massive stars. 

To date, the strongest evidence that many SHBs 
and LSBs arise from a different progenitor population 
comes from analysis of their respective host galaxy asso- 
ciations. The host galaxies of long-duration GRBs have, 
universally, been observed to have blue colours, sub-solar 
ISM metallicities, and strong emissio n features associated 
with high specific star-formation rates (ISavaglio et alj 12009] ; 
IStanek et al"]l2006l : iModiaz et alJliooi : iBergerl l2009h . The 
burst positio n, when well-constr ained, is nearly always at 
small offset (|Bloom et alJ I2002T ) an d typically traces th e 
brightest regions of the host galaxy ()Fruchter et al.l [2006l. 
which itself is typically b l ue and morphol o gically disturbed 
i|Wainwright et all I2007T I. iFruchter et ail (|2006h also find 
that LGRB host galaxies have lower luminosities on av- 
erage compared to the galaxy population probed by sur- 
veys at similar redshifts. In contrast, the host galaxies of 
short-duration GRBs to date have been observed to be 
much more heterogeneous, including both star-forming and 
non star-forming hosts. Afterglow offsets range from neg- 
ligible to many times the half-light radius of the putative 



host llProchaska et al.ll2006l ; iBloom et al1l2007l ; lBergerll2007i ; 
iTroia et al.ll2008T l. 

It is against this backdrop that GRB 090426 enters the 
scene. With an observed duration of Tgo = 1.28 ± 0.09 s 
(§2), based on its observer-frame duration alone it groups 
more closely with the SHB class, an identification that is 
further bolstered when its high redshift of z = 2.609 (JU 
is considered, implying a rest-frame duration of only 0.35 s 
that is unambiguously within the range of classical SHBs. 
By contrast, the most convin cing host association s for SHBs 
are at z < 1 (although see iBerger et alJ 120071 and a dis- 
cussion of the potentially high-redshift SHB GR B 060121; 
Ide Ugarte Postigo et ai1l2006l ; lLevan et al.ll2006l ). 

Irrespective of the phenomenological classification or 
the physical origin of this event, we report on an optical 
spectrum of the afterglow — the first ever reported for an 
event with an observed duration of < 2 s — which shows ev- 
idence of an environment quite unlike that of most (but not 
unprecedented among) GRBs of long duration with spet- 
roscopic observations. We also present the results from our 
campaign of late-time imaging and spectroscopy, which iden- 
tify the highly UV-luminous host galaxy of this event. All of 
the values and conclusions in this paper are consistent with 
our GCN Circular^] but should be considered to supersede 
our previous work on this event. 



THE DISCOVERY AND CLASSIFICATION 
OF GRB 090426 



At 12: 48:47 on 2009 April 2< 
(BAT; iBarthelmv et al " 



the Burst Alert Telescope 
onboard the NASA Swift 



satellite ([Gehrels et al.ll2004h was triggered on GR B 090426 



The X-Ray Telescope (XRT; lBurrows et al.|[2005h began ob- 
serving the field at 84.6 s afte r the trigger, and th e ultravi- 
olet/optical telescope (UVQT: iRoming et al1l2005h followed 
at 89 s after the trigger. UVOT detected a candidate optical 
afterglow at a = 12 h 3 6 m 18 3 .07, 6 = +3 2 °59 / 9".6 (J2000), 
which was reported bv lCummings et al.l (|2009T j 13.8 min af- 
ter the burst trigger. The optical counterpart at these co- 
ordinates was a lso confirmed 43.5 min after the burst by 
IXin et all (2009) based on observations obtained 76 s after 
the burst with the Teramo-Normale Telescope. The Sloan 
Digital Sky Survey (SDSS) shows no object near the position 
of the afterglow; the closest object is a faint and extended 
source at a = 12 h 36 m 19 s .49, 8 = +32°59'05".5 (J2000), 18" 
away from t he optical afterglow w ith a photometric redshift 
of 2 ~ 0.3 (jD'Avanzo et alJl2009h . As detailed in §4.1, we 
obtained a spectrum of the afterglow 1.1 hr after trigger, in- 
dependently discovering the optical afterglow by inspection 
of the guider and acquisition fr ames and determinin g an ab- 
sorption redshift of 2 = 2.609 jLevesque et al.ll2009T ). Later, 
IThone et all (|2009l ) confirmed this afterglow detection and 
redshift with a Very Large Telescope (VLT) spectrum ob- 
served 12.3 hr after the trigger. 



2 See lBloom et al] (|2008h and lZhang et al] j2007t) for a discussion 
of GRB classification, both physical and phenomenological. 



3 The GCN system [http: / /gen. gsf c .nasa.gov/ is managed and 
operated by Scott Barthelmy. 

4 UT dates are used throughout this paper. 
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2.1 The Short-Duration/Hard-Spectrum 
Bifurcation in Swift 

The question of which phenomenological classification to as- 
cribe to GRB 090426 is obviously an important one. To 
do so, we examine the hardness-duration distribution of 
Swift and BATSE (from which the original classification 
s cheme was d e rived) . Here, we measure hardness by fitting 
a iBand et alj £l993) model to the BAT spectrum and ex- 
tracting the best-fit vF v peak energy -E pe ak. Figure [1] (top 
panel) displays the durations and hardnesses for 398 Swift 
GRBs detected by Swift between December 2004 and April 
2009. Plotted in the background are hardnesses and dura- 
tions for 1728 historic GRBs detected by the BATSE exper- 
iment, taken from the current BATSE catalogu«[f]. 

Our methodology for determining - Speak (and Tgq) from 
the BAT data is described in detail bv lButler etafl J2007h. 
The p rocedure involves prior assumptions on the Band et al.l 
(| 19931 ) model parameters in order to overcome challenges as- 
sociated with fitting this broad-band spectral model to data 
covering only the BAT 15-150 keV spectral band. Large 
error bars on -E pea k necessarily result, as shown in the fig- 
ure. To estimate i5 pca k for the largest possible number of 
BATSE GRBs, we determine a relationship between the 
BATSE ca talogue hardness rat io (HR) and the measured 
Speak from iKaneko et all l|2006l ) for 325 GRBs in common. 
We find £ pea k « 80(HR 3 4i2) a69 keV, with a scatter of 0.3 
dex. Here, HR3412 is the hardness ratio of fluences in BATSE 
bands 3 + 4 over 1 + 2. 

It has been not ed previously (e.g., ICurtis et al.ll2006l ; 
IZhang fc Cho]|2008l ) that there is only weak evidence for a 
distin ct short-dura tion class in the Swift sample considered 
alone. iBandl (|2006l ) suggests that the discrepancy arises pri- 
marily as a result of the Swift increased sensitivity (relative 
to BATSE) to long GRBs, which tends to make detected 
short-duration GRBs a factor ~ 3 less common in Swift rel- 
ative to BATSE. We explore this possibility in detail here 
by correcting the observed Swift number distributions for 
sensitivity. This is accomplished using the sensitivity curves 
as a function of -E pea k and T90 duration from Figures 3 and 
4 of Band (2006). The cur ves we use a s sume an exponen- 
tial burst light curve and a IBand et all (|l993h model with 
a = — 1 and (3 = —2. 

We first fit a double elliptical Gaussian model to the 
observed distributions in T90 and -B pea k from BATSE. As- 
suming a Gaussian shape entails making the fewest assump- 
tions on the true underlying distributions, because a Gaus- 
sian is the maximum-entropy distribution in the case of 
known mean and variance (e.g.. lGregorvll2005l ). We employ 
a Markov Chain Monte Ca rlo algorithm based o n the data 
augmentation algorithm in Ivan Dvk et all |2001) to propa- 
gate errors and marginalise over the thirteen Gaussian pa- 
rameters defining the best-fit, two-class model shown in blue 
in Figure [1] (top panel). We begin by stochastically dividing 
the BATSE observations between classes, given an initial 
guess for the Gaussian parameters and also samples for the 
T90 and -Epcak values from their respective best-fit distribu- 
tions (assumed Gaussian). With this division in place, we 
find the best-fit Gaussian parameters again and draw sam- 
ples for each from the posterior distribution using the Gibbs 



sampling technique (e.g. Ivan Dvk et al.ll200ll ). The process 
is iterated, allowing us to determine 10 3 samples for each 
parameter after dropping 100 samples ("burn in") to allow 
the chain to converge. 

Next, we scale the best-fit double-class model by the 
relative sensitivity curve to obtain the contours in black, 
which are those expected for Swift. To do the scaling, we 
must assume a relation for the number of bursts gained (lost) 
as the sensitivity is decreased (increased). We assume that 
the number scales as the relative sensitivity to some power r\. 
This rj is the slope of the cumulative number d ensity versus 
flux r elation (i.e., the log iV-log S relation, as in lPreece et al.l 
2000). We take 77 = —0.75. Our resulting black curves do 
appear to better match the Swift Tgo and -E pca k distributions 
(Figure [1] top panel). 

To quantitatively judge the validity of the corrected 
model, we determine the rate increase/decrease factor for 
each GRB as a function of -E pC ak and T90 in the Swift sam- 
ple and generate from these a corrected T90 histogram. The 
Kolmogorov-Smirnov (KS) test probability that the uncor- 
rected distributions from Swift and BATSE are the same 
is 10 -6 ; however, with application of the relative sensi- 
tivity function to adjust the rate for each event, the KS- 
test probability is only 2.7%. Hence, the distributions be- 
come only marginally inconsistent with no ad hoc changes. 
More precise tweaking, which would utilise the exact spec- 
tral/temporal properties and a detailed simulation of the 
Swift trigger algorithm, would likely improve the consis- 
tency, although this analysis would be very challenging to 
conduct and is beyond the scope of this work. 



2.2 Classification by Tgo and i5 pca k for 
GRB 090426 

We download the raw, unfiltered Swift BAT data for 
GRB 090426 from the Swift Archivs0. Our reduction of these 
data to science-quality light curve s and spectra using stan- 
dard Swift tools are detailed by iButler et"ail l|2007h . We 
employ calibration files from the 2008-10-26 BAT release. 
The BAT signal in the 15-350 keV band consists of a sin- 
gle, narrow emission spike of duration T90 = 1.28 ± 0.09 s 
(T50 = 0.48 + 0.06 s). The spectrum in the 15-150 keV band 
is well modelled (x 2 jv = 54.18/45) as a single power law 



with photon index (3 = — 2.02_ 



and an energy fluence 



(15-350 keV) of 2.5^1 x 10~ 7 erg cm -2 . Using our Bayesian 
methodology ijButler et al.ll2007f ) to extrapolate to an ap- 
proximately bolometric energy release in the 1-10 4 keV band 
(source frame), we find E- mo = 4.2lojj x 10 51 erg, with a vF v 
spectral peak energy of £ pca k = 451^ keV (observer frame). 

Above, we map a 2-class model (Gaussian G\ for the 
short/hard class and Gaussian G2 for the long/soft class) 
justified based on comparing BATSE data to the Swift sam- 
ple. We will now apply this classification to Swift GRBs and 
to GRB 090426 in particular. Importantly, the precise factor 
dictating the relative Swift/BATSE number distribution at 
a given value of -E pea k and Tgo — which we derive approxi- 
mately above as arising solely from variations in the satellite 
sensitivities — does not enter into this calculation. We only 
need to know the ratio of the probabilities, which means 



5 http:/ /www. batse.msfc.nasa.gov/batse/grb/catalog/current 6 ftp:/ /legacy.gsfc. nasa.gov/swift/data. 
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the factor drops out of the relative classification calculation. 
In principle, the relative rate factor could also depend on 
variations with redshift of the intrinsic source populations 
(ignored above) at fixed values E pea k and Tgo, but we make 
the simplifying assumption here that this can be ignored. 

It is important to demonstrate that our derived Ep ea k 
and Tgo values from Swift are sufficiently similar to those 
derived from BATSE. Because we have defined E pe ak in a 
similar fashion for both experiments, this then primarily be- 
comes an issue of comparing Tgo values derived in different 
bandpasses, whereas we know GRB spectral evolution tends 
to make a given event appe ar longer when measu red in a 
lower energy bandpass (e.g.. iFenimore et al.l[l995h . Fortu- 
nately, the Swift (15-350 keV) and BATSE (50-300 keV) 
bandpasses are similar, and we can directly measure any bi- 
ases in Swift T90 values calculated in the 15-350 keV band 
as opposed to the 50-300 keV band. Considering 411 Swift 
GRBs, we find that the median decrease in Tgo when consid- 
ering the 50-300 keV band instead of the 15-350 keV band 
is only 3.8%. The decrease is < 30% for 90% of the sample, 
which is typically (> 75% of the time) contained within our 
la error bar on Tgo. Therefore, we expect systematic varia- 
tions in Tgo with energy band to not affect our classification. 

Using the Markov Chain derived for the model division 
as a function of Tgo and E P eak, we can directly determine 
the probabilistic class association for Swift GRBs taking into 
account errors in Tgo and E pca k. Using a Markov Chain for 
this purpose effectively treats the error bars on all quantities 
and allows us to marginalise over the parameters describing 
the 2-class model. The marginalisation is important, because 
there is strong overlap in the observed BATSE distributions 
which translates to uncertainty in the Gaussian model pa- 
rameters defining our BATSE classification. To classify the 
Swift GRBs, we sample 10 3 values f or Tgo and E pea k from the 
distributions in lButler et al.l (|2007T ). Each of these samples is 
used to evaluate one of the Gi/(Gi + G2) draws above. The 
class probability, which is the Bayesian evidence in favour 
of the membership in the SHB class as compared to the 
LSB class, is calculated as the median of the Gi/(Gi + G2) 
samples. 

The ratio of Gaussians Gi/(Gi + G2), evaluated at a 
particular value for Tgo and E pca k, defines the probability 
that a given burst will belong to the SHB class under our as- 
sumptions. In the bottom panel of Figure [TJ the solid black 
curve shows the projection onto the abscissa of the solid 
distribution in the top panel. This curve is the classification 
marginalised over E pca k- Note that because we have assumed 
a sensitivity correction as a function of Tgo and -Epcak at each 
value of Tgo and -B pe ak, the correction drops out of the ratio 
Gi/(Gi +G 2 ), and we can apply the BATSE Gi/(Gi +G 2 ) 
model to Swift without needing to account for relative sen- 
sitivity. Dashed curves are also shown to display the E pea k 
dependence of the curve at two -Epcak values (E pca k = 1 
MeV and E pca k = 20 keV, respectively). 

In black circles in Figure [l] (bottom panel) , we display 
the probability for each Swift burst plotted in the top panel. 
We also plot as red circles the host-frame Tgo z = Tgo /(1 + z) 
values for 138 GRBs with measured spectroscopic redshifts. 
Considering the range of observed E pea k values, we find that 
a GRB is short/hard at > 90% confidence if Tgo < 2.2 s, 
or Tgoz < 0.8 s. These limits can be used in future stud- 
ies to select Swift GRBs belonging with high confidence to 



the short/hard class. Note that our spectroscopic redshift of 
z — 2.609 for GRB 090426 enables us to derive a rest-frame 
duration Tgo z of 0.35 ± 0.03 sec. 

We find the probability that GRB 090426, highlighted 
and circled in yellow in Figure [T] belongs to the short/hard 
class is 92.8%. Even so, we must stress that such a high con- 
fidence indication could occur by chance given a sufficient 
number of detected SHBs. For ~ 400 Swift LSBs detected 
to date, the chance probability of detecting one or more 
long/soft GRBs with durations short enough and/or hard- 
ness high enough to appear short/hard with such high con- 
fidence in our scheme is > 90%. This marks a fundamental 
shortcoming in the classification by high-energy properties 
alone, where the parameter distributions suffer broad over- 
lap, motivating further investigation into the afterglow and 
host properties. It is possible that addition al high energy 
indic ators (e.g., a "lag" consistent with zero. lUkwatta et al.l 
2009) may be useful for classification, but we do not inves- 
tigate these here. 



3 ENERGETICS AND AFTERGLOW 

As mentioned above, from our Bayesian model of the 
burst parameters, we calculate an isotropic energy release 
of 4.2 x 10 erg. Until very recently, this value of Ei so 
would be considered exceptionally large for a SHB. Indeed, 
a low Ejsc is naturally expected from most SHB models 
|Panaitescu et aill2001f ). which are typically assumed to col- 
limate their ejecta less effi ciently than the collapsar model 
l|Bergerll2007l : lNakarll2007r i . However, the recent GRB 090510 
(spatially associated with an e mission-lin e galaxy) had a 
very large E iso = 3.8 x 10 52 erg (Rau;2009), suggesting that 
short-duration bursts are indeed capable of arising from very 
energetic (and/or tightly collimated) explosions as well, and 
are probably visible substantially beyond z = 1, even if they 
are not expected to be common. 

Combining the (limited) set of observations of the after- 
glow of GRB 090426 from published circular^, we find the 
optical light curve is well described by a single power-law 
decay with index ao ~ 0.8 from t < 100s until i>4x 10 4 s. 
This is similar to the inferred X-ray decay index, ax » 0.9, 
from the online compilation of N.R.B[f|. Combined with the 
derived optical to X-ray spectral index for this interval, 
Pox ~ 0.9, and assuming s tandard synchrotron afterglow 
theory fe.g.. lSari et al.ll"l998r h these results suggest an after- 
glow with a shallow electron index (p « 1.8) and a cooling 
frequency u c below the optical bandpass. If this is indeed 
the case, we can use the limit on the cooling frequency to 
constrain the parsec-scale circumburst density. Assuming a 
constant density medium, th e cooling frequency f alling be- 
low the optical requires (e.g.. iGranot fc Sarill20o3 ) : 

n>0.05 £ - 3/2 £-^ 5 2 2 , (1) 

where n is the circumburst density (cm -3 ), es is the frac- 
tion of the shock energy partitioned to the magnetic field, 
and -Eke, 52 is the kinetic energy of the outgoing blastwave 



7 See http://gcn.gsfc.nasa.gov/gcn3_archive.html 

8 http://astro.berkeley.edu/~nat/swift; see iButler et al ] l |2007h 
for details. 



© 2009 RAS, MNRAS OOO.rflfTTl 



GRB 090426 5 



(10 52 erg). Given a maximal e_e of 1/3 at equipartition, and 



10 erg, we derive a lower limit of 



these lines by fitting each line with a Gaussian using splot 
in 1RAF. 

We find that the ionised absorption lines in our spec- 
trum are saturated, which limits us to determining conser- 
vative lower limi ts for the column densities of these lines 
jProchaskal 120061 ) based on the relation between EW and 
column density for saturated lines l)Cowie fc Songaila 1986). 
We generally find lower limits for all the saturated columns 
on the order of 10 14 cm -2 . However, we are able to calcu- 
late an upper limit for JVhi based on the absence of strong 
damping wings in the Ly-a absorption feature. From fit- 
ting the line with a Voigt profile, we find an upper limit of 
Perlev et al.ll2009h or the long GRB 050911 N m < 3.2 x 10 19 cm" 2 . Our values for EW and the various 



with -Eke ;$ 

n > 0.1 cm -3 . A similar result can be derived for the case o f 
a wind-like circumburst medium fe.g.. lChevalier fc Lif l999): 
A, > 0.01 (where p = 5 x lO 11 ^* gem , chosen to corre- 
spond to a mass loss rate of M = 10~ 5 MQyr" 1 and a 
wind speed of v w = 1000 kms' 1 ). While we caution that 
this result is based on a relatively sparsely sampled op- 
tical light curve, it is clear that presence of a relatively 
bright and slowly fading optical and X-ray afterglow distin- 
guish GRB 090426 from the extremely low-density circum- 
burst environm ents inferred for th e short GRB 080503 (n < 
5 x lQ- 6 cm~ 3 ;[P. 
jPage et al.ll2006f > 



4 OPTICAL AFTERGLOW 
SPECTROPHOTOMETRY 

4.1 Observations and Reductions 

We obtained an optical spectrum of the afterglow of GRB 
090426 usi ng the Keck Low -Resolution Imaging Spectrome- 
ter (LRIS; lOke et al.lll995l ) at 13:55 on 2009 April 26, ~1.1 
hr after the BAT trigger. The observations were conducted 
in photometric conditions. We obtained two 300 s exposures 
on the LRIS blue side using the long 1" slit mask, the 680 
dichroic, and the 300/5000 grism. We observed internal flat- 
field lamps as well as spectra of Hg, Ne, Ar, Cd, and Zn 
comparison lamps to be used for wavelength calibration. 
We also obtained a 60 s spectrum of the spectrophotometric 
standard HZ 43. The observations of the GRB 090426 after- 
glow were conducted at a high airmass of 3.05; HZ 43 was 
observed at an airmass of 3.60. 

The data were reduced using IRAF0. We used the 
lrisbias IRAF task distributed by the W. M. Keck Ob- 
servatory to subtract overscan from the LRIS images, and 
apply a wavelength correction based on our internal lamp 
observations. The spectrum was extracted using an opti- 
mal extraction algorithm, with deviant pixels identified and 
rejected based upon the assumption of a smoothly varying 
profile. We flux calibrated the data using observations of HZ 
43 to derive a sensitivity curve which was then applied to the 
GRB 090426 afterglow observation. Finally, we corrected for 
a heliocentric velocity of —16.88 km s _1 and corrected the 
spectrum to rest-frame wavelengths. Our spectrum is shown 
in Figure [2] 



4.2 Analysis and Interpretation 

In our analysis of the afterglow spectrum, we initially ob- 
served a set of absorption features at 4387 A, 5030 A, 
5061 A, and 5592 A. We identify these features as Ly-a, 
Si IV A1394, Si IV A1403, and the blended C IV AA1548, 
1551 doublet at a common redshift of z = 2.609. At this 
redshift we are also able to identify the N V AA1239, 1243 
doublet, Si II A1260, and C II A1334 absorption features. 
We determine the rest-frame equivalent widths (EWs) for 



9 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under cooperative agreement with the NSF. 



column densities are given in Table [T] 

The value of the neutral hydrogen column is very low 
in comparison to other GRBs: based on the cumulative 
distribution of Nm in 28 long-duration GRBs at z ^ 2 
|Chen et al.l I2007T ). we find that the afterglow of GRB 
090426 has a lower N m than ~90% of GRB afterglow spec- 
tra. Our GRB 090426 afterglow spectrum also appears to 
have weaker low-ionisation absorption (Si II, C II) than 
~95% of previous afterglow spectra. This sets GRB 090426 
apart as atypical when compared to the host environments of 
other G RBs, which generally ha ve much stronger absorption 
features ijProchaska et al1l2008l ). Nevertheless, even among 
"typical" long-duration GRBs, such very low columns are 
not completely without precedent, and a few long-duration 
GRB afterglow spectra are found to have similarly low Nh i 
to GRB 090426. Typically, GRB afterglows with Ly-a lines 
have column densities of Nn i w 10 21 cm -2 ; one notable ex- 
ception is GRB 021004, which has N H i « 1 x 10 19 cm" 2 . 
It is suggested that the low measured Nn i in that after- 
glow spectrum is due to ionisation of t he H I by the radia - 
tion field of the massive-star progenitor l|Fvnbo et alll2005l ). 
Another unusual afterglow is associated with the long/soft 
GRB 060607; with A H i = 6.3 x 10 16 cm" 2 , it has the lowest 
H I column density of any GRB afterglow. The GRB 060607 
spectrum lacks any detection of the N V lines, though it does 
show C IV and Si IV a bsorption at the redshift of the GRB 
(jProchaska et al1l2008h . However, no host galaxy has been 
detected for GRB 060607 thus far, d own to an ff-band lim- 
iting magnitude of AB(H) = 26.5 (|Chen et all 120091 '). By 
contrast, we do in fact detect the N V doublet in the after- 
glow spectrum of GRB 090426. N V is thought to originate 
in the immediate circumburst environment of the GRB, and 
this absorption feature is quite typical of mo st other ob- 
served GRB afterglows (jProchaska et al]|2008h . 

Similarly, examples of systems with extremely weak 
low-ionization lines, while quite rare, are not unprecedented 
among ordinary long bursts: GRB 070125 and GRB 071003 
were both fou nd to have extremely weak host Mg I I absorp- 
tion systems l|Cenko et al.ll2008l : [Perlev et al.ll200Sft . indica- 
tive of a low-density galactic environment, possibly in a tidal 
tail or halo. 



5 THE HOST GALAXY OF GRB 090426 
5.1 Imaging 

On 2009 May 21 we imaged the field using GMOS-S on 
Gemini-South and the i-band filter for 20 exposures of 180 s 
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each, for 1 hr of total integration time. Images were pro- 
cessed using archival twilight flatfields and fringe corrected 
within the gemini package in IRAF. The following night 
(2009 May 22) we acquired additional imaging in the V band 
using the FOCAS instrument on Subaru. A total of 9 images 
of 300 s each were acquired for a total integration time of 
45 min. Images were processed using standard techniques in 
IRAF. Both optical images show a bright, extended object 
with complicated morphology (a bright, elongated object 
with fainter lobes of emission to the NE and S) near the 
afterglow location (Figure [3J. 

Finally, on 2009 May 31 we imaged the field using NIRC 
on the Keck I telescope. A total of 31 exposures of 1 min (10 
coadds x 6 s) were acquired in the K band, plus 9 in the H 
band (also 10 x 6 s), and 9 in the J band (3 x 20 s). Images 
were processed and stacked using a modified Python/ pyraf 
script originally written by D. Kaplan. No object consistent 
with the optical band is detected in any filter. Based on a 
calibration to 2MASS standards observed in frames taken 
later in the night, we place 3er limiting magnitudes on the 
host-galaxy flux of J > 23.0, H > 22.1, and K > 22.0 mag 
(Vega). 

To calculate the offset of the afterglow relative to the 
putative host galaxy, we aligned both the LRIS acquisi- 
tion image (taken the night of the burst) and the Subaru 
V-band observation to reference stars in the Sloan Digital 
Sky Survey, giving a position of a = 12 h 36 m 18 s .052, S = 
+32°59'09.14" (J2000). This position places the afterglow 
within 0.1" (800 pc in projection) of the centre of the north- 
eastern lobe of the system that we subsequently identify as 
the host galaxy complex. 

Aperture photometry of the brightest (central) region 
of the host as well as the knot at the afterglow location 
was performed with IRAF using a 1" radius. The resulting 
photometry, corrected for the modest Gal actic extinction 
(E(B-V) =0.017 mag; ISchlegel et al"1ll998h , is presented in 
Table [5] In addition to the spatial coincidence, the identical 
colours strongly suggest that the two objects are physically 
relate d. Interpolating to the flux at 1700 A (see lReddv et al.l 
2008), the photometric magnitude of the northeast compo- 
nent of the host corresponds to a rest-frame UV luminosity 
of approximately 0.7 L*, or ~ 2L* for the entire host com- 
plex, indicating a luminous host galaxy. 

5.2 Spectroscopy 

We obtained an additional late-time spectrum at the after- 
glow location with LRIS on 17 June 2009. Extrapolating the 
early-time optical light curve, the afterglow flux should have 
faded sufficiently such that any emission would be domi- 
nated by host-galaxy light. Our observations consisted of 
two 1500 s exposures. The blue side was configured with the 
600/4000 grism, providing coverage of 3500-5500 A with a 
scale of 0.62" pixel -1 , and a spectral resolution of ~ 4 A. 
The red side employed the 400/8500 grating with wavelength 
coverage of 5500-10,000 A, a scale of 1.18" pixel -1 , and a 
spectral resolution of ~ 7 A; however, we do not discuss 
the red-side spectrum here, since we are interested primar- 
ily in a detection of Ly-a emission from the host galaxy. The 
long, l"-wide slit was oriented with a position angle of 41.3° 
to capture both the compact "knot" at the afterglow loca- 
tion and the nearby extended galaxy, while an atmospheric 



dispersion corrector was ut ilised to account for differential 
refraction (Firi ppenkdll982l ). We are confident that the slit 
was at the correct location, because we detected another 
object at its expected position along the spatial axis of the 
spectrogram. 

The spectra were reduced in a manner identical to that 
described in §3.1. We find no sign of any flux, either con- 
tinuum or narrow emission lines, at the location of the af- 
terglow or host complex. At z — 2.609, Ly-a would fall at 
Aob a = 4389 A. Using observations of the stand ard star Feige 
34 (|Massev et al.lll988l ; I5"ke1ll990l ; IStonelll996l ) from earlier 
in the night, we place a limit on any emission- line flux at this 
location of F < 7 x 10 -17 erg cm -2 s -1 (assuming the line 
was narrow enough to be unresolved in our spectra). Us- 
ing th e sta r-formation rate conversions from iBrocklehurstl 
l|l97ll ) and iKennicuttl |l983j), we therefore place a limit on 
the unobscured star-formation rate (SFR) at the location of 
the afterglow and host of SFRly-q < 4Mq yr - . This value 
is significantly less than that derived from the k-corrected 
rest-frame UV (1500A) continuum emission (neglecting ex- 
tinction correcti ons), where SFRuv = 14.4 ± 2.0 Mq yr -1 ; 
l|Kennicuttlll998h . However, we note that SFRs derived from 
Ly-a emission can often underestimate the true SFR by 
over an order o f magnitude due to r esonant scattering, 
dust absorption ()Mas-Hesse et all 12003! '). and a strong de- 
pendency on the a ge of the star- forming population (e.g., 
IValls-Gabaudlll993h . Furthermore, the night of these spec- 
troscopic observations was not photometric (variable, thin 
cloud cover), and therefore our flux calibration may be in 
error (though likely at less than the 50% level). 



5.3 Models and Interpretation 

We generated synthetic photometry in our measured filters 
at z = 2.609 using the irregul ar, Sc/d, Sb/c, and ellip tical- 
galaxy templates in hyper-Z llBolzonella et all 12000J; tem- 
plates originally from IColeman et al. 1980h . The templates 
were screened by varying amounts of host-galaxy dust, both 
with and without the 2175 A bump (assuming an LMC 
and SMC extinction law, respectively) to compare with the 
observed colours. The elliptical and Sb/c models were im- 
mediately ruled out as incompatible with the blue i — K 
colour implied by the NIRC non-detections, as were large 
amounts of dust extinction in any case. The highly starburst- 
dominated irregular template (plus a small amount of ex- 
tinction, Av ~ 0.4 mag) is favoured over the slightly more 
evolved Sc/d template, though given the limited photome- 
try available and the simplified nature of the modelling this 
conclusion is less robust. From this examination, it is clear 
that the broad-band photometry indicates a stellar popula- 
tion dominated by young stars. 



5.4 Associating the Galaxy Complex with GRB 
090426 

In the optical afterglow spectrum we find no detections of 
any intervening absorption systems at a redshift of z < 
2.609. Similarly, in our spectroscopic observations of the pu- 
tative host complex at the position of the afterglow we see 
no spectral features that would be consistent with contribu- 
tions from a foreground system. 
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In the absence of an emission-line redshift of the com- 
plex, it is reasonable to ask what the possibility is that 
this host association is the result of a chance alignment be- 
tween GRB 090426 and a foreground system at z < 2.609 
(at ^complex much larger than z a b s the system would be 
too bright intrinsically). Examining only the northeastern 
lobe of the host co mplex and following the prescription in 
iBloom et all (|2002h . we estimate a probability of chance 
alignment between afterglow and the central region of its 
host galaxy (using an effective radius of 0.25") to be 0.1%. 
More conseratively, if we instead consider the entire host 
complex (approximately 1.8" radius), the probability of 
chance alignment is still low at 4%. Based on the low like- 
lihood of a chance alignment and a lack of spectroscopic 
evidence supporting the presence of a foreground system, 
we conclude that this is indeed very likely the host galaxy 
(complex) of GRB 090426. 



6 CONCLUSIONS 

The small astrometric offset from what appears to be a blue 
host galaxy initially seems to be difficult to square with the 
inference of a very low column density implied by the ab- 
sorption spectrum. However, it is noteworthy that the upper 
limit probes only the neutral gas; the implied UV luminos- 
ity from the bright host system suggests a large ionising 
radiation field in and around the Galactic disks which may 
have ionised a significant fraction of the neutral gas along 
the line of sight to the GRB. Furthermore, both the detec- 
tion of N V and the significant circumburst density implied 
by the bright afterglow indicate that the immediate envi- 
ronment of GRB 090426 is not dramatically different from 
those of long-duration GRBs in general. 

On the other hand, a genuine halo environment seems 
unlikely. Some degenerate merger scenarios involve a signifi- 
cant (> 1 Gyr) delay between initial formation of the system 
and the merger, suggesting that the positional and temporal 
coincidence of the afterglow with what could be a starburst 
induced by tidal interaction with the nearby object would be 
relatively unlikely. In these scenarios, the progenitor system 
is also subject to a systemic velocity "kick" during binary 
evolution that results in sign ificant linear motion of the sys- 
tem away from its birt hsite (iFrver et al.ll 19991 : IBloom et al.l 
1 19991 : although see also lBelczvnski et alj2002f7 . For instance, 
a binary with a 100 km s _1 kick perpendicular to our line 
of sight that takes 1 Gyr to merge will travel 100 kpc from 
its birthsite; in our adopted cosmolog\F°l at z = 2.609 this 
amounts to an angular distance of 12.8", compared to the 
observed < 0.2" offset. 

This does not immediately reject the association of this 
event with merger-product progenitors; if such progenitors 
can merge over a range of timescales (including relatively 
short ones), the association of SHBs with active starbursts 
would be no surprise. If GRB 090426 is interpreted as aris- 
ing from a merger, this event may suggest that SHBs may 
very well be akin to Type la supernovae (which appear to 
be generated by bo th long and short production channels; 
ISullivan et alB ooe ) . Indeed, many SHBs to date have shown 
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little to no ap preciable offset fr om their (sometimes blue) 
host galaxies (jTroia et al.l 120081 ) . This event also serves as 
a spectroscopic example of the high-redshift short-d uration 
GRB population inferred from spatial associationsin IrBergerl 
(2007). While the most direct evidence for a degenerate 
merger remains the detection of a gravity wave signature 
(see for example [Bloom et alj|200gh . then GRB 090 426, at 
a large red shift with large Ej ao , would suggest (cf. iBergerl 
l200l 2009) unfortunately that a significant number of SHBs 
detected by BATSE and Swift occur well outside of the Ad- 
vanced LIGO volume. 



The simplest conclusion from the available observations 
of the afterglow and host galaxy is that GRB 090426 is 
more closely linked with the core collapse of a massive star. 
The implications of this association are no less profound: 
they indicate that the mechanism that generates gamma 
rays is capable of operating on timescales as short as 0.3 s, 
imposing strong demands of the central engine; in the most 
basic collapsar model for GRBs, the duration timescale is 
generally as sumed to be at least an order of magnitude 
longer (see IWooslev fc Blooml 120061 ) . While events like 
GRB 090426 are probably rare (due to relative volumetric 
effects), these inferences also cast significant doubt on 
the classification of a large population of what would 
otherwise have been considered classical short-duration 
bursts: if this burst had occurred at a similar redshift to 
prototypical SHBs 050509B or 050724 (at z = 0.2-0.3) it 
would have fallen unambiguously within the SHB duration 
distribution. This also illustrates the insufficiency of T90 
alone as a classification criterion, given the 92.8% likeli- 
hood that GRB 090426 is a member of the short/hard 
phenomenological class. At minimum, we feel that at low 
redshift, the search for accompanying supernova emission 

- an unambiguous sign of a genuine massive-star origin 

— remains vital to properly distinguishing among different 
progenitor scenarios. 
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Figure 1. Top: The T90 durations and -E pca k values for 398 Swift GRBs. In the background, we plot (small blue circles) 1728 T90 and 
.Epeak values for historic GRBs detected by the BATSE experiment. Overplotted is the best-fit double-Gaussian model to the BATSE 
data (blue; 50% peak probability and 5% peak probability contours). The black curves show the relative distortion expected for these 
distributions appearing in Swift, given the relative satellite detection efficiencies (see Band 2006). Because of its increased relative capacity 
to trigger on long/soft events, there are relatively fewer (by a factor of ~3) short/hard events compared to long/soft events in the Swift 
sample as compared to the BATSE sample, making classification given the Swift data alone difficult if not impossible. GRB 090426 is 
highlighted and circled in yellow. Bottom: The ratio of Gaussians defines the probability that a given burst will belong to the short/hard 
class. Red circles give the host-frame Tgo z = Tgo/(l + z) values for 138 GRBs with measured spectroscopic redshifts. The solid black 
curve shows the projection onto the abscissa of the solid distribution in the top panel. Dashed and dotted curves are also shown to 
display the £ pea k dependence of the curve at two £p ea k values (i? pca k = 1 MeV and -E pca k = 20 keV, respectively). GRB 090426 Tgo 
and Tqo z are circled in yellow. 
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Table 1. Species Detected in the Keck/LRIS GRB 090426 
Afterglow Spectrum 



Species (A ) EW (A) Nx (cm" 2 ) 



Ly a (1215.67 A) 


2.8 


± 


0.1 


< 


3.2 


X 


10 19 


N V (1238.82 A) 


0.7 


± 


0.1 


> 


2.8 


X 


10 14 


N V (1242.80 A) 


0.3 


± 


0.1 


> 


1.8 


X 


10 14 


Si II (1260.42 A) 


0.6 


± 


0.1 


> 


3.8 


X 


10 13 


C II (1334.53 A) 


0.3 


± 


0.1 


> 


1.0 


X 


10 14 


Si IV (1393.75 A) 


2.2 


± 


0.1 


> 


3.2 


X 


10 14 


Si IV (1402.77 A) 


1.7 


± 


0.1 


> 


3.7 


X 


10 14 


C IV (1548.20 A/1550.78 A) 


3.6 


± 


0.1 


> 


9.1 


X 


10 14 



Note. — EWo and Ao are given in rest-frame quantities. 



Table 2. Photometry of the GRB 090426 Host-Galaxy 
Complex 



Filter 


Date 
(2009 UT) 


Telescope /Instrument 


Extended Host 
(AB Magnitude) 


Compact Knot 
(AB Magnitude) 


V 


May 22.26 


Subaru / FOCAS 


24.21 ±0.15 


24.73 ±0.15 


a 


May 21.05 


Gemini South / CMOS 


24.09 ±0.15 


24.61 ± 0.18 


j 


May 31.30 


Keck I / NIRC 


> 23.9 


> 23.9 


H 


May 31.30 


Keck I / NIRC 


> 23.5 


> 23.5 


K 3 


May 31.35 


Keck I / NIRC 


> 23.8 


> 23.8 
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Figure 2. Keck spectrum of the GRB 090426 afterglow. The spectrum was observed with LRIS on Keck I at 13:55 on 2009 April 26, 
~1.1 hr after the BAT trigger. The observations were conducted in photometric conditions. The data were reduced using IRAF, and have 
been corrected for a heliocentric velocity of —16.88 km s -1 . We plot both the observed wavelength (lower abscissa) and the rest-frame 
wavelength at our redshift of 2.609 (upper abscissa). We note detections of the Ly-ct, N V, Si II, C II, Si IV, and C IV features at this 
redshift. 
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Figure 3. False-colour optical image of the host-galaxy field from combined i-band data from GMOS-S on Gemini South and V-band 
data from FOCAS on Subaru. A magnified region of the host complex is inset at top right. The afterglow position identified by our LRIS 
acquisition imaging is shown in both images as a yellow circle of radius 0.2" (2a ) and is consistent with the northeast component of the 
complex. The large galaxy 18" to the East of the host complex is that noted bv lDAvanzo et al.l ||2009|) . 
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